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Abstract. The Poynting-Robertson (P-R) effect has been applied to meteoroids for sev¬ 
eral decades. It is well known that the P-R effect produces only changes in the orbital 
plane of the particle. The differential equations governing the secular changes in both 
eccentricity and semi-major axis are known since the time of Robertson, provided the 
initial orbits are not near circular. A disadvantage of this type of osculating orbital ele¬ 
ments is that it cannot be applied to motion of arbitrarily shaped particles. Relevant type 
of osculating orbital elements is discussed. It’s application to the P-R effect is presented 
and simple analytical formulae for secular changes in both eccentricity and semi-major 
axis are derived. 


1. Introduction 

The relativistic equation of motion for a perfectly absorbing spherical dust particle under 
the action of electromagnetic radiation was obtained by Robertson (1937). Robertson 
derived expressions for the secular changes in the orbital elements. The effect essentially 
reduces both semi-major axis and eccentricity and is popularly known as the Poynting- 
Robertson (P-R) effect. The results of Robertson were verified by Wyatt and Whipple 
(1950) and applied to the evolution of meteoroid streams. One particularly important 
aspect, first pointed out by Harwit (1963), is that a meteoroid escapes from the Solar 
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System when the orbital energy becomes positive and that this can happen when the 
energy due to the radial component of the radiation force is included, without it being 
necessary for this force to exceed the gravitational attraction. 

A relativistic generalization corresponding to the P-R effect was discussed by Klacka 
(1992a). A more simple derivation is given in Klacka (1993a). There have been many 
misunderstanding, and even errors, in the treatment of radiation force effects in the lit- 
terature and some of these are discussed in Klacka (1993b, 2000b). We will not dwell 
further on most of these here. A consistent set of equations governing, to the hrst order 
in u/c, the secular changes in the orbital elements were given by Klacka (1992b) (here, v 
is velocity of the particle, c is the speed of light), where the appropriate initial conditions 
were also discussed. This work, and most of the preceeding work assumes that the orbits 
are not near circular. When the orbits are near circular, different considerations must ap¬ 
ply, for example a circular orbit can not reduce in radius without increasing in eccentricity, 
contrary to the normal behaviour under P-R effect. The near-circular situation has been 
discussed by Klacka and Kaufmannova (1992, 1993); the results analytically conhrmed 
Breiter and Jackson (1998) and their nonphysical conclusions are discussed in Klacka 
(2000a, 2001a). We will not consider this near-circular case in detail, but discussion will 
contain a short comment. 

The Poynting-Robertson effect also causes an advancement of perihelion, although 
the the effect is second order in v/c (that is, the secular change is proportional to c“^). 
Correct result can be found in Balek and Klacka (2002). 

The aim of this paper is to discuss relevant form(s) of osculating orbital elements and 
their secular changes. The motivation comes from the fact that the P-R effect represents 
only a very special case of interaction between a particle and incident electromagnetic 
radiation (Klacka 2000c, 2000d, 2000e, 2001b, 2002a, 2002c, Klacka and Kocifaj 2001). 
Type of secular changes of osculating orbital elements derived by Robertson (1937) cannot 
be compared with secular changes of these elements for real particles: secular changes of 
osculating orbital elements based on central gravity force alone, are relevant physical 
quantities (Klacka 2002b). The correct equations governing the secular evolution of semi¬ 
major axis and eccentricity are derived for the P-R effect. 


2. P-R effect: Equation of motion 

Consider a spherical particle with a geometric cross section A' in a flow of incident 
radiation with flux density S' in a direction denoted by a unit vector S. Let Q),,, denote 
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the efficiency factor for radiation. As is usual within a relativistic context, we use 7 = 
1/a/1 — !(?. It is convenient to define w = 7(1 — v ■ Sjc). 

The relativistically covariant equation of motion for the P-R effect may be expressed 
as (Klacka 1992a, 2000c, 2002a, 2002c) 


d T 




( 1 ) 


where is four-vector of the particle of mass m, p^ = m u^, four-vector of the 
world-velocity of the particle is = {"f v) and = (l/w) (1, S). 

Relativistically covariant form of equation of motion is important in deriving an ex¬ 
pression for the secular change of an advancement of perihelion, which can be found in 
Balek and Klacka (2002). Moreover, higher orders in v/c play an important role when 
making calculations for near circular orbits (compare Breiter and Jackson 1998, Klacka 
2000a, 2001a). This will not be considered here, except section 3.1. Thus, reduction of 
Eq. (1) to the first order in v/c will be used, if it is not explicitly stressed. 

Applying Eq. (1) to the situation where a particle is moving in gravitational and 
electromagnetic fields of a star (usually the Sun) with mass M gives 



( 2 ) 


where p, = GM, G being the gravitational constant, (3 is the ratio of radiation pressure 
force to gravity force, that is /? = A' r'^ S / [c p m). 


3. Secular changes of orbital elements radiation pressure as a part of central 
acceleration 


We have to use — ^ (1 — P) S / as a central acceleration determining osculating 
orbital elements if we want to take a time average (T is time interval between passages 
through two following pericenters) in an analytical way 



assuming non-pseudo-circular orbits and the fact that orbital elements exhibit only small 
changes during the time interval T; ap is semi-major axis, ep is eccentricity, fjs is true 
anomaly, pp = ap{l — e|); the second and the third Kepler’s laws were used: dfp/dt = 
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a//z(1 — f3)pi3 - conservation of angular momentum, = ^(1 — /3)/(47r^). (For more 

details see Klacka, 1992b.) 

Rewriting Eq. (2) into the form 

vS 

c 



iif 

at ( 


we can immediately write for components of perturbation acceleration to Keplerian 
motion: 


^0 R = - 2 p — - , Fp T = - P ^ - , Fj3 N = 0 , 

C C 


( 5 ) 


where Fj^ r, Ffj x and F )3 ^ are radial, transversal and normal components of pertur¬ 
bation acceleration, and two-body problem yields 


vr = 


Vt = 


f3) 


P0 


P0 


e/3 sin fa , 


(1 -f epcosfp) . 


( 6 ) 


The important fact that pertubation acceleration is proportional to v/c (<C 1) ensures 
the above metioned small changes of orbital elements during the time interval T. 

Perturbation equations of celestial mechanics yield for osculating orbital elements 
{ap - semi-major axis; ep - eccentricity; ip - inclination (of the orbital plane to the 
reference frame); Vlp - longitude of the ascending node; top - longitude of pericenter; Qp 
- 0^ = top + fp is the position angle of the particle on the orbit, when measured from 
the ascending node in the direction of the particle’s motion): 


dap 

dt 


2 ap 
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1 - ej \l - (3) 
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cos 1(3 


( 7 ) 


where r = pp/{l + ep cos fp). 
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Inserting Eqs. (5) - (6) into Eq. (7), one easily obtains 

da /3 _ ^ M 2 a /3 1 + e| + 2 e /3 cos ffi + e| sin^ 

dt c 1 — e| ’ 

^ (2e/3+ e/3sin^//3 + 2 cos//3) , 

f ' 

^ 

dt 

duj/3 II 1 1 f \ ■ f 

—JT = -(2 - COS//3 Sin/^ , 

at H c e /3 

dQf3 ^ \/m(1 - P) Pl3 

dt ~ r2 


( 8 ) 


It is worth mentioning that da/j/dt < 0 for any time t (we again stress that near circular 
orbits are not considered). 

The set of differential equations Eqs. ( 8 ) has to be complemented with initial condi¬ 
tions. If the subscript 0 denotes orbital elements of the parent body (e. g., comet), then 
even if a particle is ejected with zero relative velocity, it will move on a different orbit 
due to the central force on it being reduced because of the radiation. (For more general 
case, with non-zero ejection, see Gajdosik and Klacka 1999.) The elements of the ’’new” 
orbit are given by 


0/3 in = ao (1 - /3) 1 - 2/3 


, 1 -I- eo cos /o 

1 - 
i Cq 


-1 


_ /i 1 - eg - 2/3 (1 -f eo cos /o) 

^/3 in — \ ^ 


a-py 

where /o = ©o — wq, W/j has to be obtained from 


(9) 

( 10 ) 


e/3 in COS (©0 ai/3 in) 
e/3 in sin (©0 0)0 in) 

f2/3 in — , '10 in — 


_ /3 -f eo cos /o 
1-/3 

_ eosin/o 
1-/3 ’ 

fo 7 ©/3 in — ©0 ■ 


( 11 ) 

( 12 ) 


Some figures may be found in Klacka (1992b). 

By inserting Eqs. ( 8 ) into Eq. (3), taking into account that e /3 /„ < 1 (see Eq. (10)), 
one can easily obtain the secular changes of orbital elements: 


(^) = -/3^ 

at c 


2 -f 3e^ 


/de0 


dt ' 2 ^ c 


0/3 


(1 - ej) 


3/2 


e/3 


o? 


- 4 ) 


1/2 ’ 


(13) 


(14) 
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( 


d&f3 

dt 


\/m (1-/3) 

3/2 


5.i. Advancement of perihelion, v'^/-terms 


(15) 

(16) 


Working to first order, as above, in v/c yields no secular change of to (advancement of 
perihelion in Solar System). Thus, we have to include higher orders. Including now terms 
up to je?, Eq. (1) yields 



which gives 


^ ji VRVT 

R 2c 2 ^r2 c2 


(17) 


(18) 


In reality, these components of perturbation acceleration represent only a part of the 
total perturbation acceleration, since general relativity theory produces other terms. The 
total advancement of pericenter is given (compare Balek and Klacka 2002) as 

^^^^■c2af (l-e2) (1-/3)^/^ ■ 

It can be easily verified that < dujp / dt > is i) an increasing function of (3, ii) the 
perihelion circulates in a positive direction, iii) the rate of the advancement of perihelion 
is not bounded for /3 ^ 1. 


4. Electromagnetic radiation and general equation of motion 


General equation of motion of a particle under the action of electromagnetic radiation 
can be written in the form (Klacka 2002c) 


dp^^ 


= E 

1=1 


wlS A' 


O' 


cl) W 


/ c) 


( 20 ) 


where also nonradial components (in the proper frame of reference of the particle) of 
radiation pressure are considered and also thermal emission of dust particle is taken into 
account {dm/dr = 0). 


Let us consider orbital evolution of real dust particle, under the action of gravitational 
and electromagnetic radiation fields of the central body (star). We can write for the Solar 
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System (Klacka 2002c) 


d V 

d t 
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_ d 2 1*00 
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G M© m c Jo 

fd2 = ^ ® / i?0(A) { - GL,(A/«;) ffi(A/rc)} JA , 

G M© m c Jo 

^3 = { - ^-a(A/w;) ff^(A/u;)} JA , 

G Mq me Jq 

w = 1 — v-ei/c. 


( 21 ) 


( 22 ) 


Rq denotes the radius of the Sun and Bq{X) is the solar radiance at a wavelength of A; 
G, Mq, and r are the gravitational constant, the mass of the Sun, and the distance of the 
particle from the center of the Sun, respectively. The asymmetry parameter vector g' is 
defined by g' = (l/G^^^j) / n'{dC'g^g^/dXl')d^', where n' is a unit vector in the direction 
of scattering; g' = g[ e[ + g'^ e'^ + e'^, e[ = {1 + v ■ ei/c) ei - v/c, e' • e' = Sij. 

The first set of Eqs. (22) is important from the point of view that general equation of 
motion represented by Eqs. (20) and (21) has to be reduced to the special case treated 
by Einstein (1905). 

We see that the motion of real, arbitrarily shaped, particle may by much more com¬ 
plicated than that corresponding to the P-R effect; P-R effect is a special case of Eq. 
(20). To the first order in v/c, Eq. (21) reduces to Eq. (2) for (32 = /d^ = Q/i = Q'e 2 = 
Qg3 = 0; (di = (d, when comparing Eq. (2) and Eq. (21). Motion of a real particle depends 
on optical properties of the particle, e. g., shape, chemical composition of the particle, 
mass distribution within the particle (porosity). 


5. P-R effect and secular changes of orbital elements gravitation as a central 
acceleration 

If one would like to inspire with calculation od orbital elements as it was presented in 
section 3, he should consider — /i(l — (di) S / as a central acceleration determining 
osculating orbital elements even for general case represented by Eq. (21). However, the 
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quantity (3i depends on optical properties of the particle and, thus, also on particle’s 
position with respect to the source of electromagnetic radiation. 

It is not wise to use — ^ (1 — /3i) S' / as a central acceleration determining oscu¬ 
lating orbital elements for general case represented by Eq. (21), since ^ (1 — (ii) changes 
almost randomly during a motion. Thus, it is not wise to use — /i (1 — j3) S / as 
a central acceleration determining osculating orbital elements for the P-R effect, if we 
want to compare the evolution of orbital elements for the P-R effect and the evolution of 
orbital elements for general case represented by Eq. (21). We need something which is not 
changing almost randomly. The wise quantity is /i. Thus, central acceleration — /i S / 
will be used as a central acceleration determining osculating orbital elements. 


5.1. P-R effect and perturbation aquations of celestial mechanics 

We use — S / as a central acceleration determining osculating orbital elements. 

Taking into account Eq. (2), we take the action of electromagnetic radiation as a 
pertubation to the two-body problem. We can immediately write for components of per¬ 
turbation acceleration: 


Ffl = /34-2/34 — > Ft = -P 


ii il 

r ’2 /-> 


Fn = 0 


(23) 


where Fr, Ft and Fn are radial, transversal and normal components of perturbation 
acceleration, and two-body problem yields 


VR 


Vt 


I — esinf , 

p 


— (1 -f ecosf) . 
P 


(24) 


Perturbation equations of celestial mechanics yield for osculating orbital elements (a 
- semi-major axis; e - eccentricity; i — inclination (of the orbital plane to the reference 
frame); ft - longitude of the ascending node; co - longitude of pericenter; Q - Q = to f 
is the position angle of the particle on the orbit, when measured from the ascending node 
in the direction of the particle’s motion): 


da 

dt 
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dt 

di 

dt 

dQ 

dt 
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1 — e^ 


{Fr e sin / -b Ft (1 -b e cos/)} 
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^/JFp 


Fr sin / -b Ft 
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e -b COS / 
1 -b e cos / 


sin 0 
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du> 


1 


P 


— = - - J- {FrcosJ-Ft 


e 

r 


dt 


y/WP 


P 

^ sin© 

Pn —.—r COS I , 


2 + e cos / 
1 + e cos / 


sinf} - 


sinz 

r 


sin 0 

Fn —r COS I 


^/JPp sin i 


where r = p/{l + ecosf). 

Inserting Eqs. (23) - (24) into Eq. (25), one easily obtains 
da 


11 

„2 


pL a'- 


dt 
de 
dt 
di 
dt 
dQ 
dt 
duj 

dt r 

de _ y/JTp 

dt ~ r 2 


. . ^^ 2 al + e^ + 2 e cos / + sin^ / 

2 e sm/ - j3— - 


1 - e 2 .•' "'r 2 c 1 - e 2 

= /3 sin / — - ( 2 e + e sin^ / + 2 cos /) , 

j,Z ^ \ / 

= 0 , 

= 0 , 

= - (i i cos/ - /3 (2 -ecosf) sin / , 

" 6 C 6 


(25) 


(26) 


It is worth mentioning that da/dt < 0 for any time t does not hold (perturbation corre¬ 
sponds to complete nongravitational acceleration, and, thus, Eqs. ( 8 ) do not hold). 

The set of differential equations Eqs. (26) has to be complemented with initial condi¬ 
tions. If the subscript 0 denotes orbital elements of the parent body (e. g., comet), then 
particle ejected with zero relative velocity will have initial osculating orbital elements 
identical with those of the parent body: 


— CIq , Cin — Cq 5 — ^0 j 

^in — flo , tOin — ^0 5 ©m — ©0 ■ (2"^) 

The important fact is that Eqs. (26) contain also terms not proportional to v/c (<C 
1). These important terms protect us to use procedure analogous to that represented 
by Eq. (3). While dispersion of osculating orbital elements is very small during a time 
interval T for the case when pL [1 — (3) is used in central acceleration, the dispersion 
of osculating orbital elements may be large during the same time interval for the case 
when pi is used in central acceleration. Thus, any formal averaging of Eqs. (26) leading 
to equations analogous to Eqs. (13)-(16) is not correct. This fact was discussed in Klacka 
(1992b), see also Figs. 1 and 2 in Klacka (1994). 


5.2. Expressions for secular changes of semi-major axis and eccentricity 

We have already explained that it is not allowed to make a simple time averaging anal¬ 
ogous to that described by Eq. (3), when — pi S / r"^ is used as a central acceleration 
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determining osculating orbital elements. However, it is of interest if we have to numer¬ 
ically solve Newtonian vectorial equation of motion (Eq. (2)) and make numerical time 
averaging (over a time interval between passages through two following pericenters), or 
some analytical simplifications can be done - something analogous to Eqs. (13)-(14). 

Correct answer is: yes, we can make some analytical simplifications when we want 
to obtain secular changes of semi-major axis and eccentricity even when — fj, S / is 
used as a central acceleration determining osculating orbital elements. We will derive the 
correct equations in the following two subsections. 


5.2.1. Radial forces and orbital elements 

We will proceed according to Klacka (1994), in this subsection. 
Let us consider a gravitational system of two bodies 


V = -j , 


(28) 


where en = ei = S. Let perturbation acceleration exists in the form 


a = f3^eR 


0 < P < 1. Thus, the final equation of motion is 


V = — 


{1 - P) 


br . 


(29) 


(30) 


Eq. (30) yields as a solution the well-known Keplerian motion and the orbit is given 


by 


Pc 


1 -I- ec cos (0 — ojc) ’ 
where 


(31) 


Pc = Oc (l - el) . 


(32) 


The subscript “c” denotes that orbital elements are constants of motion. If we write 


V = vr er -|- vt er , 


(33) 


where bt is a unit vector transverse to the radial vector br in the plane of the 
trajectory (positive in the direction of motion), we have 


VR = \J p - P) Pc^ Be sin (0 - LOc) 


(34) 


VT = y P ip - P) Pc^ [I + Be COS (0 - UJc)] ■ 


(35) 
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In principle, we may consider also a new set of orbital elements, which are defined by 
the central gravitational acceleration. Equations (31), (32), (34) and (35) are then of the 


form 


_ P 

1 + e cos (0 — w) ’ 

(36) 

p = a (l - e^) , 

(37) 

vr = 0 ^ p~^ e sin (0 — w) , 

(38) 

vt = VP P~^ [1 + e cos (0 - w)] ; 

(39) 


the fact that 0 is unchanged in both sets of orbital elements is used. 

Position vector and velocity vector define a state of the body at any time. Equations 
(31) and (36) yield then 


Pc 


(40) 


1 + Cc COS (0 — ujc) 1 + e cos (Q — uj) ' 

Analogously, the other two pairs of equations (Eqs. (34) and (38), and, Eqs. (35) and 
(39)) give 


\J{I - P) Pc ^ Cc sin(0 - uJc) 
^{I - (}) pc^ [1 + Cc cos(0 



e sin (0 



- e cos(0 


w)] . 


(41) 

(42) 


One can easily obtain, using Eqs. (40) and (42), 


Pc {I - P)=p , 

and Eqs. (41)-(42) yield then 


(43) 


(1 — (3) Cc sin (0 — LOc) = e sin (0 — w) , (44) 

(1 — /3) [1 + Cc cos (0 — ojc)] = 1 + e cos (0 — w) . (45) 

Eq. (45) can be written in the form 

(1 — (3) Cc cos (0 — ojc) — j3 = e cos (0 — w) . (46) 

Eqs. (44) and (46) yield 

= (1 — (3)^ + 0 ^ — 2 j3 {1 — (3) Cc cos(0 — ujc) ■ (47) 


Equation for u is given by Eqs. (44) and (46), using also Eq. (47). Finally, Eqs. (32), 
(37), (43) and (47) yield 
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a = Qc <1 + P 


1 + + 2 Cc cos (0 — LUc) 

1 — e'i 


(48) 


Eqs. (47)-(48) show that orbital osculating elments a and e change in time, during an 
orbital revolution - the larger /3, the larger change of a and e. 

The osculating orbital elements e and a obtain values between their maxima and 
minima, which can be easily found from Eqs. (47)-(48). One can easily verify that these 
relations hold: 


emin = 1(1 - P) Sc - p\ <e< {1 - P) Cc + p = e, 

\ — Ba .CL . 1 + e 


a. 


min 


^ a ^ i H- a-n 

ac I - Cc + P (I + ec) ~ ac ~ I + ec + P (I - Cc) ac 


(49) 

(50) 


5.2.2. Mean values of semi-major axis and eccentricity 

Eqs. (49) and (50) represent interval of values for eccentricity and semi-major axis, when 
— fj, S / is used as a central acceleration determining osculating orbital elements. How¬ 
ever, we can make time averaging during a period T, which was decribed by Eq. (3): 


(e) = 

1 

1 


T Jo ’ 

(a) = y 

1 ^ 
P (1 

^ _ 

d - P) Pc 

dt 

> j^2 

4 


a{t) dt 


Pc 


1 -I- Cc cos/c 


Eqs. (51) yield 
(e) = (1 - e^) 


2\3/2 1 


e(/c) 


2 TT Jo (1 + Cc cosfcY 


= (1 - 


^277 


Yfc) 


2 tt Jo (1 + Cc cosfcf 
Using Eqs. (47) and (48), we finally obtain 


dfc , 


dfc- 


3/2 1 a/(1 - Pf el+P"^- 2P (1 - P) ec cos /c 

(e) = (1 - e,) — - - - df, , 

Jo [1 + Cc cos fc) 

3/2 J_ [l+ p{l + el + 2ec cos fPj / (l - eg)] 

27r 


z) = ac(l-e^)'- 


/o (1-I-CcCOS/c) 

The following properties can be verified: 
i) (e) > /?, (e) = P for Cc = 0; ii) (e)/ec > 1, (e) = Cc for /3 = 0 
iii) 9(e)/9ec > 0; iv) d{e)/dp > 0; 

v) (a) > ac/ (1 -I- P), (a) = ac/ (1 -I- P) for Cc = 0; 

vi) {a)/ac < 1, (a) = Oc for P = 0; 

vii) d{a)/dec > 0; viii) d{a)/dp < 0; ix) d{a)/dac > 0. 


dfc 


(51) 


(52) 


(53) 
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5.2.3. Secular changes of semi-major axis and eccentricity 


Summarizing our results, it is possible to calculate secular evolution of eccentricity and 
semi-major axis, according to the following prescription. 

At first, initial conditions for 0^3 and 6,3 are calculated: 

, 1 -f eo cos /o 

1 - 

1 Cq 


(a/3)m = ao (1-/3) 1-2/3 




1 - - 2/3(1 -I- eo cos/o) 


(54) 


{i-py 

supposing that particle was ejected with zero ejection velocity from a parent body - 
quantities with subscript 0 belongs to the parent body’s trajectory. 

As the second step, the set of the following differential equations must be solved for 
the above presented initial conditions: 


^ ^ _ o _ 

dt c 

ap 

de/3 _ 5 

dt 2 c 


2 -h 3e^ 


(i-ej) 


3/2 ’ 


e/3 


-4) 


1/2 • 


(55) 


Finally, semi-major axis and eccentricity are calculated from: 


= a/3(l-e^) 


3/2 1 


/3 ^1 -I- -I- 2 e /3 cos a;^ / ^1 - 


1 -1 


27r 


e = (1 - e|) 


3/2 1 

27r 


3 (1-I- 6/3 cos x) 

-^(1 - /3)^e| -I- /32 - 2/3(1 - /3) e /3 cosx 

(1 -|- 6/3 COSx)^ 


dx 


dx . 


(56) 


The set of equations represented by Eqs. (54)-(56) fully corresponds to detailed numer¬ 
ical calculations of vectorial equation of motion, if we are interested in secular evolution 
of eccentricity and semi-major axis (supposing (e/ 3 )i„ < 1 and 6/3 does not correspond to 
pseudo-circular orbit) for the case when central acceleration is defined by gravity alone. 

It is worth mentioning that instantaneous time derivatives of semi-major axis and 
eccentricity may be both positive and negative (see Eqs. (26)), while secular evolution 
yields that semi-major axis and eccentricity are decreasing functions of time. 


6. Discussion 

When the orbits are near circular (pseudo-circular) - central acceleration contains ra¬ 
diation pressure term - the orbit can not reduce in semi-major axis without increasing 
in eccentricity. Both types of orbital elements, defined by central accelerations, have 
been used in detail numerical calculations in papers by Klacka and Kaufmannova (1992, 
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1993). Due to the property (e) > (3 (see section 5.2.2), one must be aware that even 
value (e) > (3 may correspond to pseudo-circular orbit. The results for pseudo-circular 
orbits were analytically confirmed by Breiter and Jackson (1998). Although the analyti¬ 
cal approach reproduces known results without detail numerical calculations, it paralelly 
produces nonphysical results which may not be distinguishable from the correct results. 
The nonphysical analytical results are caused by usage of P-R effect in the first order in 
V jc - very special analytical solutions will diminish when higher orders in v jc are used. 
The nonphysical analytical results have been discussed in more detail by Klacka (2000a, 
2001 a). 

7. Conclusion 

We have discussed orbital evolution for the P-R effect. General equation of motion for 
interaction between particle and incident electromagnetic radiation shows that radiation 
cannot be considered as a part of central acceleration. The central acceleration has to 
contain only gravity of the central body (star). In order to compare secular changes of 
semi-major axis and eccentricity for real particle and the P-R effect, the paper presents 
secular changes of these orbital elements for the P-R effect (Eqs. (54)-(56)). 
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